Gamma Ray Astronomy with Magnetized Zevatrons 
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Nearby sources of cosmic rays up to a Z6V(= 10^^ eV) could be observed with a multi-messenger 
approach including secondary 7— rays and neutrinos. If cosmic rays above ~ 10^* eV are produced in 
magnetized environments such as galaxy clusters, the flux of secondary 7— rays below ~ 1 TeV can 
be enhanced up to several orders of magnitudes compared to unmagnetized sources. A particular 
source of enhancement are synchrotron and cascade photons from e^e~ pairs produced by protons 
from sources with relatively steep injection spectra oc E~'^'^. Such sources should be visible at the 
same time in ultra-high energy cosmic ray experiments and 7— ray telescopes. 
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Introduction. One of the central unresolved issues 
of modern astroparticle physics is the origin of cosmic 
rays, particularly those at Ultra High Energy (UHECR) 
which have been observed at energies up to a few times 
10^° eV Q. Sources capable of accelerating such parti- 
cles, like powerful radio galaxies commonly found inside 
galaxy clusters and groups, are thought to be rare and 
have yet to be identified. 

Astroparticle physics is currently experimentally 
driven and involves many different existing or planned 
projects ranging from UHECR observatories such as the 
Pierre Auger project to neutrino telescopes 0|i well 
as ground and space based 7— ray detectors operating at 
TeV and GeV energies, respectively 0. It is clear that 
GeV-TeV 7— ray and neutrino astronomy will prove an 
invaluable tool to unveil the sources, and probe into the 
mechanism, of UHECRs. Even if a putative source were 
to produce exclusively UHECRs, photo-pion ^ and pair 
production by protons on the cosmic microwave back- 
ground (CMB) would lead to a guaranteed secondary 
photon and neutrino fluxes that could be detectable. 

Secondary photon fluxes from UHECR interactions 
with the CMB have been discussed before in the lit- 
erature 0j Hi 0- these works, however, photo-pair 
production by protons, photons generated by the GZK 
interaction and the effect of structured magnetic fields 
on both the electromagnetic (EM) cascades and UHE- 
CRs were not taken into account in a quantitative and 
comprehensive way. In Ref. |lCll| proton acceleration up 
to ^ 10^^ eV around cluster accretion shocks was stud- 
ied which, however, can not explain the highest energy 
cosmic rays. In the present letter we find that in combi- 
nation, these ingredients can increase secondary photon 
fluxes from individual UHECR sources below ~ 1 TeV up 
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to several orders of magnitude. This potentially makes 
their detection much easier. 

Furthermore, we extend the investigation to also con- 
sider the case of steep proton spectra oc E~^-^ at ultra- 
high energies. This is motivated by the scenario where 
extragalactic protons dominate the observed flux down 
to the "second knee" at ~ 4 x 10^^ eV, such that the an- 
kle at ~ 5 X 10^^ eV is caused byphotopair production 
by the extragalactic protons [lllfl^ . 

In the following we compute the expected 7-ray flux 
from 10 MeV to the highest energies, by combining 
simulations of UHECR propagation in structured large 
scale extragalactic magnetic fields [3 (EGMF) with 
production of secondary hadrons by nucleon primaries 
using the event generator SOPHIA jl^]. Pair produc- 
tion by protons on the CMB is taken into account as a 
continuous energy loss: a proton with energy E gener- 
ates electron-positron pairs (heretoafter simply referred 
to as electrons) with a power-law energy distribution 
dn/dEe oc E~''^'^ (ov Ee < E All the electromag- 

netic products of these interactions are then followed to 
the observer using an EM cascade code based on Ref. 0| , 
which takes into account the inhomogeneous distribution 
of magnetic fields in the simulation box and the presence 
of a cosmic infrared background from 17]. 

We consider the case of a discrete source in one of the 
prominent magnetized galaxy clusters from the simula- 
tions based on Refs. J^, with a size comparable to 
the VIRGO cluster. Sources in such clusters could be ac- 
tive galaxies (AGNs) whose hot spots, for example, have 
been suggested to accelerate UHECRs 0. We use c = 1 
throughout the letter. 

Influence of Pair Production and Proton Injection 
Spectra on GeV-TeV Photon Fluxes. We assume point 
sources whose flux contributes a fraction 77 < 1 to 
the total solid integrated UHECR flux observed around 
IQi^eV, i.e. ~ 2.5 x 10^i9(ry/0.01) cm'^ s'^, the ap- 
proximate average over existing flux measurements [l|. 
For distances d < 500 Mpc this flux is not too much in- 
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fluenced by energy losses and roughly correponds to an 
UHECR injection power above 10^^ eV of 
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FIG. 1: Comparison of differential 7-ray fluxes (multiplied by 
squared energy) generated by GZK interactions and photo- 
pair production for a source at 100 Mpc distance injecting 
protons up to 1 ZeV with a spectrum oc E^"^'"^ (solid and 
dashed-dotted lines) and with a spectrum oc E~^'^ (dotted 
line) in the absence of magnetic fields. The power emit- 
ted above 10^^ eV is Lig = 10*"^ erg s~^, corresponding to 
77 ~ 0.02 in Eq. Also shown are point flux sensitivities 
of the 7-ray instruments EGRET [23, GLAST [13, HESS 
2 [23, and MAGIC The dashed fine on the right is the 
observable primary UHECR flux from the source. 



Fig. n demonstrates the contribution of pair produc- 
tion to the photon fluxes at GeV and TeV energies. It 
increases with the steepness of the proton injection spec- 
trum because pair production is the dominant energy loss 
process for protons with energies 10^* ^ ^ 4 x 10^^ eV. 
It appears that, for steep enough injection spectrum 
oc E^"' with a ^ 2.6 — 2.7, necessary if extragalactic 
cosmic rays dominate above the ankle |llL Il2l | , and large 
enough distances between the source and the observer, 
photo-pair production dominates over pion production 
for the secondary 7-ray flux. 

We now investigate the consequences of the low energy 
extension of the UHECR injection spectrum. 

Cosmic ray protons of energy E and integral flux 
Jq^{E) confined within a volume ~ which interact 
with a baryon gas of density rib produce 7— rays of en- 
ergy ~ fjE ~ O.l-E at a rate ~ Jl.^{E)app nt ■ Here, 
the proton-proton cross section app ~ 3 x 10~^^ cm^ can 
be approximated as energy independent. For proton con- 
finement times iconf ^ R, we can express the confined 
integral flux in terms of the total injection rate /"^ (E) as 



J^j, (£:) ~ IciiE) tcon{{E)/R^- Furthermore, in a steady 
state situation, the leaking cosmic ray flux observed at 
a distance d is Jg^"(S) ~ M{E) P^I{E)/{4tt(P), where 
M{E) < 1 is a modification factor accounting for inter- 
action losses during propagation to the observer. We can 
then relate the integral photon fiux at energy above E^ 
from pp interactions within the volume R^, Jpp{E^), to 
the integral UHECR fiux above Eqk, JcniEcR), 



J1?{E^) ^ a-ppnbtconi{E.y/ f.y) 



(2) 

Both JPP and Jcr(^cr) are fiuxes observed at distance 
d. 

For a galaxy cluster Tif, ~ 10^'^ cm^'^, R ~ 2 Mpc and 
tconfiE-y / f-y) < 10^" yr (the age of the Universe). Thus, 
the optical depth for pp interaction is cTppUhtcoui ^ 0.3. 
Furthermore, at Ecr ^ 10^^ eV, A/(£'cr) ~ 1 and with 
the cosmic ray fiux from Eq. |^ we have 



Jf (i?,) 7.5 X 10-° (^) 



0.017 /^'j^(10i9eV) 
The upper limit on the 7-ray flux at E^ 
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100 MeV 

from EGRET is typically 4 x 10"^ cm"^ s"^ for clus- 
ters like Coma or Virgo ^2^. Eq. Q thus implies the 
condition 
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For an unbroken power law /^p(i5) oc E-^~^ for Ej/f^ ^ 
E < 10^^ eV, this would imply the relatively strong con- 
straint 



f3 < 2.17- 0.1 log(77/0.01). 



(5) 



This constraint can be avoided if the power law cuts off 
or becomes harder at low energies such that Eq. |0J is 
satisfied. In particular, in the scenario in which extra- 
galactic protons dominate down to a few 10^^ eV, their 
injection spectrum, a ~ 2.6 [llLll2| . cannot continue be- 
low ~ ioii+iog(';/o oi)/i.6gY^ 

At the same time, for /? > 2., the total power emit- 
ted by the source in cosmic rays down to energy Scr" 

is LcR ~ {lO^^eY/E'^^)'^''^ L19. Therefore, low en- 
ergy cosmic ray fiux extensions with power law index not 
much larger than 2 also assure reasonable total cosmic 
ray powers which remain largely below the high end of 
bolometric luminosity for AGNs, L < 10'** erg s^^. 

Sources in Magnetized Galaxy Clusters. Fig. |2] demon- 
strates the influence of the cluster EGMF on the fluxes 
of secondary 7— rays. Note that all fluxes scale with rj. 
This implies that for magnetized galaxy clusters and rel- 
atively soft injection spectrum oc the TeV 7— ray 
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signal should be visible at least with HESS 2, provided 
r] > 0.02, whereas detectability by MAGIC and GLAST 
requires 77 <^ 0.3. For harder injection spectra oc E~^-^ 
these numbers are 77 ^ 0.05 and rj ^ 0.5, respectively. 
These figures are for a source at d = 20 Mpc, but depend 
only moderately on d. 
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FIG. 2; Differential 7-ray fluxes (multiplied by squared en- 
ergy) from photo-pion and pair production by UHECR in- 
jected with an spectrum by a source at 20 Mpc. 
We assume 77 ~ 0.3, corresponding to a proton luminosity 
1/19 — 7 X 10*'^ erg s~^ in Eq. (Q. Compared are different 
magnetic field and propagation models, as indicated. 

Currently the isotropy of UHECRs at lO^^^ eV imposes 
only loose bounds on 77 due to the small statistics. Up- 
coming generation experiments like Auger will constrain 
77 much better in a near future. 

Fig.|21shows that the total amount of energy going into 
7— rays is roughly 10% of the energy in protons above the 
GZK threshold and does not depend much on the intr- 
acluster magnetic field. This is because, except for very 
steep injection spectra, most of the energy going into 
photons is due to pion production whose mean free path 
above ~ 6 x 10^^ eV is already smaller than the source 
distance. In contrast, the tail below ~ 1 TeV is due to in- 
teractions of cosmic-rays below the GZK threshold whose 
energy loss length is much larger than the source distance 
such that their path length can be considerably increased 
by the intracluster magnetic field. At these energies the 
dominant emission mechanism is synchrotron radiation 
from e^e~ pairs produced by UHECR protons, as we 
will see in the following. 

Electrons of energy E in a magnetic field B emit syn- 
chrotron photons of typical energy 

^--^■«^i«"(io5^)X(rniG)°^- 

The typical energy of electrons and photons produced 
in pion production is 5 x 10^^ eV 6], whereas most of 



the pair production occurs for proton energies E between 
~ 10^^ eV and ~ 4 x 10^^ eV, and gives rise to an elec- 
tron energy distribution dn/dEg oc EZ'^^'^ . Therefore, in 
a 0.1/xG field, the synchrotron emission from the elec- 
trons produced in the first stages of an EM cascade initi- 
ated by pion production occurs mainly below '^0.1 TeV. 
For pair production, if the proton injection spectra are 
steeper than E~'^ , most of the EM energy is produced 
by protons of a few times 10^* eV which end up in syn- 
chrotron photons below ~ 1 TeV with a long tail to lower 
energies due to the rather flat pair spectrum. Both these 
effects are seen in Fig. |21 

The TeV photon flux in Figs. ^ and |21 can be approxi- 
mated as 

Jf (ii;^ - 1 TeV) - 6.3 X 10-" (^) cm-^s-^ . (7) 

Requiring the TeV 7— ray flux to be dominated by the 
UHECR interactions rather than low energy interac- 
tions, JPP(~ TeV) < TeV), amounts to the con- 
dition /3 % 1.95. As a consequence, in order for both the 
EGRET constraint and a TeV 7— ray flux to be domi- 
nated by UHECR interactions rather than low energy 
interactions, would require a hard cosmic ray injection 
spectrum below ultra-high energies. 

In Fig.[21the one-dimensional simulation neglected pro- 
ton deflection and used the 3d profile of the magnetic field 
projected onto the line of sight. The 3d structure of the 
EGMF tends to enhance the photon flux from pair pro- 
duction considerably. This is because UHECR between 
the pair production threshold at ^ 10^^ eV and the pion 
production threshold can diffuse transverse to the line of 
sight producing photons at ~ GeV energies for timescales 
much longer than the rectilinear propagation time. For a 
nearby source this would result in a substantial enhance- 
ment of this photon flux. This effect could also lead to a 
GeV- TeV 7— ray halo whose structure could be observ- 
able in the case of a nearby powerful source with steep 
UHECR injection spectrum. In Fig. O the cumulative 
flux dflE^dN/dEdrt is represented for a source in a 
magnetized cluster at 20 Mpc from the observer. In this 
case the 7— ray halo has a spatial extension of order 3° 
and is dominated by pair production at angles ^ 3°. Such 
a source located at 100 Mpc (like Coma) would have an 
extension ~ 0.6°, still resolvable by imaging atmospheric 
Cerenkov detectors. 

GeV-PeV cosmic-ray protons and TeV electrons ac- 
celerated at cluster shocks in galaxy clusters can also 
produce diffuse 7— ray emission at a comparable level 
through pp and inverse Compton emissions respec- 
tively [24]. The radiation spectrum produced by these 
mechanisms is a flat power law, E'^dN/dE oc E~°' with 
a ~ 0. Thus it should be distinguishable from the spec- 
tra illustrated in Fig. 12 characterized by a broken power 
law with a < at TeV energies. Notice that the latter is 
rather insensitive to the slope of the injected UHECRs as 
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over the size of the magnetized region embedding the 
UHECR source, and the TeV source could therefore be 
spatially resolved. 
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FIG. 3: Spatial extension of the counterpart in 7— ray above 
a TeV of the magnetized source at 20 Mpc in the 3d model of 
Fig.m The relative contributions of pair production and GZK 
photons are shown separately as cumulative fluxes emitted at 
off-sets from the source center larger than 9. 



the emitting particles are produced in a cascade process. 

Finally, the photon fluxes do not depend significantly 
on iSjnax in our scenarios, provided E'max <^few IO^°eV. 

Conclusions. Ultra-high energy cosmic rays produce 
secondary 7— ray from pion production and pair produc- 
tion on the cosmic microwave and other low energy pho- 
ton backgrounds. If a significant fraction of highest en- 
ergy cosmic rays is produced in galaxy clusters which are 
known to contain magnetic fields of fractions of a micro- 
Gauss over Mpc length scales, the secondary 7— ray below 
~ I TeV could be detectable by 7— ray experiments such 
as HESS 2, and potentially also by MAGIC and GLAST. 
This is especially the case for relatively steep injection 
spectra oc E^^-^ above 10^^ eV which are required by sce- 
narios explaining the ankle by pair production of extra- 
galactic protons. Injection spectra steeper than oc E~'^''^, 
however, cannot continue down to ~GeV energies. In- 
stead, they have to become harder somewhere between 
GeV and ultra-high energies to avoid over-production of 
photons produced in inelastic pp collisions. 

Whereas the 7— ray flux around a TeV turns out to be 
relatively insensitive to the magnetic field, its slope below 
~ TeV contains information about the cluster magnetic 
field and the ultra-high energy cosmic ray injection spec- 
trum. Furthermore the 7— ray flux is expected to extend 
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